Ab initio study of the intrinsic exchange bias at the SrRuOs/SrMnOa interface 
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In a recent publication (S. Dong et at, Phys. Rev. Lett. 103, 127201 (2009)), two (related) mech- 
anisms were proposed to understand the intrinsic exchange bias present in oxides heterostructures 
involving G-type antiferromagnetic perovskites. The first mechanism is driven by the Dzyaloshinskii- 
Moriya interaction, which is a spin-orbit coupling effect. The second is induced by the ferro- 
electric polarization, and it is only active in heterostructures involving multiferroics. Using the 
SrRuOa/SrMnOs superlattice as a model system, density- functional calculations are here performed 
to verify the two proposals. This proof-of-principle calculation provides convincing evidence that 
qualitatively supports both proposals. 

PACS numbers: 75.70.Cn, 75.30.Et, 75.80.+q 
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I. INTRODUCTION 

The exchange bias (EB) effect is known to be present 
in a wide range of magnetic composites. This effect has 
been extensively used in a variety of magnetic storage 
and sensor devices In principle, the physical origin 
of the EB is the magnetic coupling at the interface be- 
tween antiferromagnetic (AFM) and ferromagnetic (FM) 
(or ferrimagnetic) materials. However, a more funda- 
mental microscopic theoretical understanding of the sub- 
tle EB effect has not been developed in spite of the fact 
that the EB phenomenon has been known for more than 
half a century and used in industry for decades. In many 
textbooks of magnetism, the EB is described as induced 
by spin pinning effects at the FM/AFM interface. An 
uncompensated AFM interface, which means there is a 
net magnetic moment at the interface on the AFM side, 
is often invoked to illustrate how the pinning may work, 
as sketched in Fig. [IJa). These uncompensated mag- 
netic moments provide a bias field to the neighboring 
FM moments via the exchange coupling across the in- 
terface. This ideal scenario is very clear with regards to 
the possible cause of the effect, but it appears incomplete 
since it actually fails in several real cases. - 

To understand the deficiencies of the uncompensated 
AFM scenario, consider the case of AFM interfaces that 
are optimally uncompensated, as in Fig. [TJa). An ex- 
ample is provided by the (111) surface of G-type AFM 
materials. In these large EB would be expected. 

However, such an effect has not been observed in real 
materials with these characteristics.^ In contrast, several 
compensated AFM interfaces, as illustrated in Figs. [IJb) 
and (c), do display an EB effect, that may be even 
larger than for the "uncompensated" interface of the 
same materialsji To solve these puzzles, several possi- 
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FIG. 1. (Color online) Sketch of possible FM-AFM interfaces. 
The upper/lower (blue/red) layers denote FM/AFM compo- 
nents, respectively. Local magnetic moments are marked as 
arrows, (a) A fully uncompensated AFM interface, as it oc- 
curs for example in an (001) surface of A-type AFM state or 
(111) surface of G-type AFM state, (b-c) Fully compensated 
AFM interfaces, as they occur in a G-type AFM state. In (a) 
and (b), the FM and AFM magnetic moments are coUinear. 
In (c), the FM and AFM magnetic moments are noncoUinear, 
namely the magnetic easy axes of the FM and AFM materials 
are different. 



ble mechanisms have been proposed in the past decades. 
Extrinsic factors are often considered, such as interface 
roughnesSj-'^ spin canting near the interface,!'^ frozen in- 
terfacial and domain pinning,- ^^'^^ and others. For more 
details, the readers are referred to the theoretical review 
by Kiwi that summarizes several previous models/^ In 
most of these models, a small "frozen" uncompensation 
of the AFM moments near the interface remains the com- 
mon ingredient for the explanation of the bias field, de- 
spite the different origins of this uncompensation. 

In recent years, considerable attention and research 
have been devoted to the study of the exchange bias 
in magnetic oxides heterostructures.—*^ In particular, 
when a multiferroic material participates as one of the 
two components at the interface, the exchange bias can 



then be tuned by electric fieldS)^"— an effect that is 
difficult to understand with only purely magnetic argu- 
ments. Motivated by these experimental studies, some of 
the present authors have recently proposed two (related) 
mechanisms to understand the EB effects in magnetic ox- 
ides heterostructures.^^ These mechanisms are based on 
the Dzyaloshinskii-Moriya (DM) interaction and on the 
ferroelectric (FE) polarization, with the later only active 
in the heterostructures that involve multiferroics. These 
two mechanisms are quite different from previous models 
used to address the EB effect, since they do not rely on 
the existence of uncompensated AFM moments anymore. 

The DM interaction was first proposed over half a cen- 
tury ago (two years after the observation of the EB effect) 
to explain the existence of weak ferromagnetism in some 
AFM materialsj^^*^ The DM interaction originates in 
the spin-orbit coupling (SOC) effect caused by relativis- 
tic corrections to superexchange. This DM interaction 
can be written as Dij ■ (Si x Sj), where Si and Sj are 
nearest-neighbor spins. The vector Dij obeys the asym- 
metric rule Dij — ~Dji. In transitional metal oxides, 
Dij is determined by the bending of the Mi-O-Mj bond 
(where Af 's denote transition metal ions and O the oxy- 
gen ion). From symmetry considerations, the direction of 
Dij should be perpendicular to the plane spanned by the 
bended Mi-O-Mj bond>2i The intensity of \Dij\ is pro- 
portional to the oxygen displacements from the middle 
point of Afi-Mj-;^ Thus, if the Mi-O-Mj bond is straight, 
namely a 180° angle bond, then Dij becomes zero. 

According to the formula for the DM interaction, the 
presence of a non-collinear spin pair is a key ingredi- 
ent. However, these spin arrangements are not very com- 
mon as ground states in bulk materials since most spin 
states are typically governed by strong superexchange in- 
teractions that often favor collinear states. However, it 
is quite natural to have noncoUinear spin pairs across the 
interfaces in heterostructures, due to the different mag- 
netic easy axes of the two components. 

In bulks of several perovskites, there are intrinsic and 
collective tiltings and rotations of the oxygen octahe- 
dra, which are often generically called the "GdFeOa- 
type" distortions;^ In several oxide thin films and het- 
erostructures these octahedral tiltings and rotations are 
still present, although the distortion modes may change 
with respect to those that exist in the bulk limit due to 
the lattice mismatch and substrate strain or stressi^"— 
A crucial aspect of our arguments is that due to these 
collective octahedral distortions, the Mi-O-Mj bonds are 
bent in a staggered manner, giving rise to the staggered 
arrangement of D vectors shown in Fig.[2](b). According 
to our recent model investigations, a uniform bias field 
emerges from the combination of the staggered D vectors 
and the staggered G-type AFM spin order 

Furthermore, if one component has an (in-plane) FE 
polarization which induces a uniform displacement be- 
tween cations and anions, the angles of the Mi-O-Mj 
bonds become staggered, as sketched in Figs, [^b) and 
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FIG. 2. (Color online) (a) Optimized superlattice crystal 
structure (drawn by VESTA) The space group is Pmc2\ 
(orthorhombic). The in-plane lattice constants a and 6 are 
fixed as 3.905v^ A to match the SrTiOa substrate. The opti- 
mized c-axis constant is 15.5 A. The pseudo-cubic x-y-z axes 
are shown. The collective rotations/tiltings of oxygen octahe- 
dra are obvious to the eye, giving rise to bended Mn-O-Mn, 
Mn-O-Ru, and Ru-O-Ru bonds, (b) The optimized Mn-O- 
Ru bond angles a are 170.4°, while the bending directions 
between nearest-neighbor pairs are staggered. The D's point 
"in" (0) and "out" (®) for these two bonds respectively, (c- 
d): When an (in-plane) FE polarization P (the interfacial 
oxygens displaced along the [110] direction by 0.005\/2 A) is 
imposed to the original bending bonds, then the bond an- 
gles become modulated: /3=168.5° and 7=171.6°. Here the 
difference between /3 and 7 is enlarged for a better visual 
comparison. 



(c). As a result, the intensities of the exchange couplings 
are modulated, which also induces a uniform bias field 
when in the presence of the G-type AFM spin order 

Therefore, these two mechanisms rely on the crystal 
distortions, the spin-orbit coupling, and the spin-lattice 
coupling, all of which are intrinsic in oxide heterostruc- 
tures. However, from the experimental information avail- 
able up to date it remains difficult to reach a definitive 
conclusion about the validity of the proposed mecha- 
nisms, particularly considering the complexity of real ox- 
ide heterostructures. Thus, it is very important to ana- 
lyze further our proposal by considering other theoretical 
perspectives to judge how robust our ideas actually are. 
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II. MODEL SYSTEM AND METHOD 

To verify the two proposals for the EB effect discussed 
above, here an approach complementary to that fol- 
lowed in the original publication Ref. [l^ will be pursued, 
namely the ab-initio density-functional theory (DFT) 
simulation will be employed. It should be remarked that 
although the DFT methods are becoming more and more 
powerful in several branches of condensed matter physics 
and materials science, it is rare to find DFT applications 
in the field of exchange bias. The reason is that in most 
previous scenarios proposed to understand the EB effect, 
large length scales are involved, and the phenomenon is 
described as induced by magnetic domains and meta- 
stable states. These subtle effects are believed to be be- 
yond the ability of current DFT techniques. However, 
the two mechanisms proposed in Ref. [12 are intrinsic and 
they are supposed to occur at the atomic scale, which 
allows to investigate their existence employing a DFT 
methodology. 

In this publication, the model system is chosen to be 
a SrRuOs/SrMnOs superlattice. This system is ideal for 
our purposes for several reasons: (1) It has been known 
from experiments that the Mn spins lie in-plane and form 
a G-type AFM order while the Ru spins point out-of- 
plane and form a FM orderj^"^ (2) Ah layers (SrO, 
Ru02, Mn02) are non-polar, which excludes the possible 
transfer of charge due to the polar catastrophe. Since it 
is expected that the Fermi energy of metallic SrRuOs lies 
within the gap band of insulating SrMnOa, then charge 
transfer due to different work functions are not expected 
to be of relevance either; (3) The A-site cations are both 
Sr^"*", which provides a unique interfacial termination. 
And the EB has indeed been observed experimentally in 
these superlattices^^^Ti^^ 

The DFT calculations were here performed based on 
the projected augmented wave (PAW)^ pseudopoten- 
tials using the Vienna ab initio Simulation Package 
(VASP).^i^ The valence states include 4s4p5s, SpAsSd, 
4p5s4d, and 2s2p for Sr, Mn, Ru, and O, respectively. 
The electron-electron interaction is described using the 
generalized gradient approximation (GGA) method. 

A superlattice structure (SrRu03)2-(SrMn03)2 (to- 
tally 40 atoms) is stacked along the [001] axis, as shown 
in Fig. mja). The initial atomic positions are set as in 
the bulk arrangement for SrRuOa and the in-plane crys- 
tal lattice constants are set as 3.905 A to match the 
widely-used SrTiOa substrate. Then, the out-of-plane 
lattice constant is relaxed. And the atomic positions of 
the superlattices are also fully optimized as the Hellman- 
Feynman forces are converged to less than 0.01 eV/A. 
This optimization and electronic self-consistency are per- 
formed under a non-magnetic state. 

After the optimization of the atomic positions, differ- 
ent magnetic profiles are applied into the system. The 
initial magnetic moments are ±3 for Mn (G-type AFM) 
and 2 for Ru (FM). The options in the codes to include 
the noncoUinear spins and the SOC will be switched on 



or off in our investigations. 

Due to its SOC origin, the DM interaction is quite 
weak, thus a high precision calculation is needed. In all 
the following calculations reported here, the PREC tag 
is set as ACCURATE. All calculations have been carried 
out using the plane- wave cutoff of 550 eV and a 9 x 9 x 3 
Monkhorst-Pack (MP) /c-point grid in combination with 
the tetrahedron method42. These high energy cutoff and 
high MP fc-grid make the present DFT calculations very 
CPU- and memory-demanding for the noncollinear spin 
states with SOC. 

Even with the accuracy described above, the signals for 
a bias field remain very weak. Inspired by the differential 
circuit for weak signals methods, the DFT energies are 
here always compared in pairs to get the difference AE, 
a procedure that can considerably eliminate systematic 
inaccuracies during the calculation. In each pair, the FM 
moments of the SrRuOs layers are flipped by 180°, while 
all other inputs remain exactly the same. The energies 
before and after the spin flipping are compared to verify 
whether a bias exists. With these methods, the system- 
atic fluctuations in our calculations can be suppressed 
below 1 meV or even lower^i*^ 



III. RESULTS AND DISCUSSION 

The optimized lattice is shown in Fig. [D^a). When 
in the bulk, SrMnOs is cubic4^ while SrRuOa is 
orthorhombic.^* Thus, the oxygen octahedral distortion 
in the SrRuOa-SrMnOs superlattice (on SrTiOa) is inher- 
ited from SrRuOs. Here aU Mn-O-Mn, Mn-O-Ru, and 
Ru-O-Ru bonds are bending: the in-plane bond angles 
are 176.3° for Mn-O-Mn and 156.2° for Ru-O-Ru; the 
out-of-plane bond angles are 173.8° for Mn-O-Mn, 167.7° 
for Ru-O-Ru, and 170.4° for Ru-O-Mn. It is clear that 
the bonds in SrMnOa are very close to straight while the 
bonds in SrRuOa are more bending, implying here that 
this bending of the bonds is due to SrRuOs, in agreement 
with the arguments presented above. Here, the tilting an- 
gle of the octahedra in SrRuOs (180 - 167.7 = 12.3°) is 
close to, but slightly different, from the previous DFT 
result (10.5°) on a pure FM SrRuOg fllm on SrTiOg-^ 
Considering the differences between the studied systems 
and the differences in the details of the calculations, this 
observed nice agreement between our present results and 
the previous ones suggests that our approach is reason- 
able. 

Let us consider first the collinear spin state without the 
SOC. The output local magnetic moments are 2.56±0.05 
/iB for Mn and 1.42 ± 0.01 /iB for Ru, which are close 
to the experimental data and to previous DFT calcu- 
lations (with LSDA or LSDA+U) in the corresponding 
bulk materials.~"^ These magnetic moments suggest the 
robust stability of the G-type AFM -I- FM spin configu- 
rations in the (SrRu03)2-(SrMn03)2 superlattice.^^ The 
weak disproportions (±0.05 for Mn and ±0.01 for Ru) 
arise from the different possible bonds (Mn(t)-O-Ru(t) 
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versus Mn(4,)-0-Ru(t), where f and I denote the direc- 
tions of the spins). 

The energy difference (AE) is only 0.058 meV, which 
is aheady far below the expected precision of VASP, since 
even the default accuracy for electronic minimization is 
0.1 meV in VASP. Thus, it is clear that these two states 
before and after the Ru's spin flipping must be degen- 
erate, or in other words no bias effect is present. This 
degeneracy is certainly to be expected considering the 
presence of independent inversion symmetries in the spin 
space and crystal space. 



A. DM driven EB 

To verify the previously proposed DM-driven EB 
mechanism, the noncollinear spins are now considered. 
The initial spins of Ru (S'ru) point out-of-plane while 
the spins of Mn (S'mii) are in-plane, as observed in the 
real heterostructures.'^^"'^^ Furthermore, the Mn spins are 
rotated within the x — y plane for comparison. 

At first, the SOC is kept "disabled" in the calculations 
(thus, no DM interaction is possible in principle). The 
Ai?'s are 0.231 and 0.264 meV when S'mh's are along the 
X and y axes, 0.0016 and —0.023 meV when the S'mii's 
are along the two diagonal directions, and both 0.169 
meV when the S'mu's are along (±2, ±1,0) and (±2,=f1,0). 
These small Ai?'s remain below the valid precision of 
VASP, implying again no bias effect. Thus, just having 
noncollinear spins is not sufficient to induce exchange 
bias. 

Note that since the global breaking condition for the 
electronic self-consistent loop is 0.1 meV in the VASP 
package, the allowed error in the total energy is also of 
this magnitude. The observed small energy differences 
comparable to or below this threshold may originate from 
the self-consistent algorithm of VASP. These differences 
are smaller than the desired precision and they will not 
affect any physical conclusion. Only those energy differ- 
ences larger than 1 meV (which is already a quite high 
precision for DPT) are considered to be physical in our 
calculations. 

In the next step in our analysis, both the noncollinear 
spins and the SOC become "enabled" in the DFT calcu- 
lation. The SOC will contribute to not only the site mag- 
netocrystalline anisotropy but also the DM interaction. 
Now the energy differences Ai5's become quite nontrivial, 
since they depend on the spins' directions, as summarized 
in Tab. M 

By keeping the S'ru spins along the z-axis, the in-plane 
rotations of Smii's give an isotropic trace: /S.E is 3.404 
meV when Sym is along the x-axis, and —4.224 meV 
when Smii is along the y-axis. The absolute values of 
these two AE^s are close since the effective precision of 
VASP is of the order of 1 meV. Ideally, these two val- 
ues should be identical considering the in-plane 90° rota- 
tional symmetry. The observed small but non-negligible 



TABLE I. DFT calculations with noncollinear spins but with 
the SOC disabled. The Smh and S-R-a vectors shown here 
are the initial magnetic moments which will be reduced to 
the self- consistent values: ~2.5 — 2.6 /xb for Mn and ~ 1.4 
/iB for Ru along their initial directions. The small canting 
components are ~ 0.2 /iB (out-of-plane) for Smu and ~ 0.1 
/^B (in-plane) for S'ru, as shown in Fig. [S^a). 





A£(meV) 


(±3,0,0) (0,0,2) 


0.231 


(0,±3,0) (0,0,2) 


0.264 


(±2,±2,0) (0,0,2) 


-0.023 


(±2,t2,0) (0,0,2) 


0.016 


(±2,±1,0) (0,0,2) 


0.169 


(±2,^1,0) (0,0,2) 


0.169 


^ > < ^ (a) 




^^ .^| 


it 






(b) 






It t 



FIG. 3. (Color online) Sketches of interfacial noncollinear spin 
configurations after the self-consistency calculation reaches 
convergence. Here the initial spins Smu's are along the (100) 
direction while Sru's are along (001). All output spins show 
some canting, but small in value, from the initial axes. The 
canting components are shown in the insert boxes, where the 
lengths of arrows are in proportional to the values of canting 
momenta, (a) Results without SOC. The canting angles of 
Smu's are about 4° while those of Sru's are about 5°. After 
the flips of Sru's (from left to right in the figure), all canting 
directions also flip symmetrically, (b) Results with SOC. The 
canting angles of Smu's are about 3° (left) or 4° (right) while 
those of Sru's are about 9° (left) or 4° (right). The canting 
angles (directions and values) are no longer symmetrical after 
the flips of Sru's, which also suggests a bias due the SOC 
effect. 



A£"s are certainly not due to the inaccuracy in the cal- 
culations, which should be below 1 meV according to the 
estimations discussed before, but they arise from a bias 
field that originates in the SOC. Even more interestingly, 
/S.E reaches a maximum (absolute) value —10.125 meV 
when Smu is along the [110]-direction but it is almost 
zero (0.339 meV) when Smii is along the [liO]-direction in 
the pseudo-cubic frame. This nontrivial behavior exactly 
agrees with the D 's vector directions, namely here the D's 
across the interfaces are along the [110] direction since the 
interfacial oxygens' (in-plane) displacement is along the 
[110] direction, as found from the optimized crystal struc- 
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TABLE II. DFT calculations with the noncoUinear spins and 
the SOC both enabled. The Syin and Sru are input magnetic 
moments. The output ones are ~2.5— 2.6 fiB for Mn and ~ 1.4 
Hb for Ru along their initial directions. For the first four cases, 
the observed small canting components are ~0.1 — 0.2 fiB 
(out-of-plane) for Smu and ~0.1 — 0.2 hb (in-plane) for 5ru, 
as shown in Fig. |3Ib). For the fifth case, the small canting 
components are < 0.1 /is (out-of-plane) for Smu and ~ 0.1 
fiB (j/-axis) for S-r.u- For the sixth case, the small canting 
components are ~0. 1 — 0.2 /ib (j/-axis) for ^Mn and --^0.1 — 0.2 
fiB (a;-axis) for Sku- 



Smu 




AE{meV) 


(±3,0,0) 


(0,0,2) 


3.404 


(0,±3,0) 


(0,0,2) 


-4.224 


(±2,±2,0) 


(0,0,2) 


-10.125 


(±2,t2,0) 


(0,0,2) 


0.339 


(±3,0,0) 


(2,0,0) 


-0.272 


(±3,0,0) 


(0,2,0) 


-0.311 



ture. Note that since SrRuOa is a FM metal with itiner- 
ant magnetic moments, the DM formula D ■ (S'ru x -S'Mn) 
which is based on localized moments may be inaccurate 
quantitatively. Even though, our DFT results confirm a 
quahtative agreement with DM ideas from the anisotropy 
of AE. 

To further confirm the presence of a DM contribution, 
the S'ru spins are also initialized to be in-plane. How- 
ever, in this case AE is almost zero no matter whether 
'S'Rull'S'Mn or S'ru -L <5Mn, as shown in the last two lines of 
Tab. ini Thus, all these DFT data provide clear qual- 
itative evidence that the bias field is in the form of 
D ■ {Sun X Sru) and D is along the [110] direction of 
the pseudo-cubic frame. 

The maximum value of the bias field driven by the DM 
interaction is found to be about 2.5 meV per Mn-O-Ru 
bond (here totally four Mn-O-Ru bonds across two in- 
terfaces due to the periodic boundary conditions), which 
is larger than our previous estimation^^ by two orders 
of magnitude. There are two reasons for finding such a 
prominent bias field. The most important one is that 
in our previous estimation a very small bond bending 
(= 1°) was assumed, which was clearly an underestima- 
tion for the present (SrRu03)2/(SrMn03)2 case (w 10°, 
as shown in Fig. (Ub)). Another possible reason for the 
discrepancy may be induced by the Ru cations because 
the SOC of the Ad electrons is much stronger than for 
the 3d ones. In any case, regardless of the actual rea- 
sons for the discrepancies in the estimations, finding a 
stronger effect than previously predicted only makes our 
proposals to understand the EB effect more robust. 



TABLE III. DFT calculations with interfacial FE polariza- 
tions. Here only coUinear spins are considered, namely the 
SOC is disabled. The Smu and Sru are input magnetic mo- 
ments. The output ones are ~2.5 — 2.6 /ib for Mn and ~ 1.4 
fj,B for Ru. 





Sb.u 


5r (A) 


AE (meV) 


±3 


2 


(0.005,0.005,0) 


16.307 


±3 


2 


(-0.005,-0.005,0) 


-16.486 



B. FE driven EB 

To fully verify our previous proposals, it is now nec- 
essary to check the bias field that is driven by the 
FE polarizations. However, here neither SrRuOa nor 
SrMnOs is a FE material. Then, to test this FE mech- 
anism it is necessary to apply an "artificial" polariza- 
tion at the interfaces. For this purpose, the oxygen 
anions between Ru and Mn are slightly displaced (Sr) 
uniformly from its optimized position to mimic an in- 
plane polarization)^ These displacements provide mod- 
ulated Ru-O-Mn bonds (7 = 171.6° and /3 = 168.5° 
when Sf = (0.005, 0.005, 0)A, while the original one is 
a = 170.4°), as shown in Figs. HJb), (c), and (d)). Al- 
though these displacements are imposed manually, their 
result, namely the presence of modulated bonds, had al- 
ready been confirmed in our previous DFT optimized 
analysis of La3/4Sri/4Mn03/BiFe03 superlattices since 
BiFeOs is a multiferroic material (for more details, see 
the supplementary material of Ref. [23 ). Since the pro- 
posed FE-driven EB does not require noncoUinear spin 
configurations, here only collinear patterns are consid- 
ered and the SOC is disabled. The results are summa- 
rized in Tab. Hill 

From the results shown above, it is obvious that the in- 
plane FE polarization at interfaces can indeed induce an 
EB field via the modulation of normal exchanges (note 
that the total effect on the DM interaction from such a 
bond angle modulation will be cancelled at the first or- 
der approximation, and thus can be neglected especially 
when the weak intensity of the DM interaction is consid- 
ered). The bias field is about 4 meV per Ru-O-Mn bond 
and the FE polarization corresponding to \6f\ = 0.005-\/2 
A is about 3.5 mC/m^. Although to be expected, it is 
very interesting to note that this EB field can be switched 
between positive and negative values by flipping the FE 
polarization using electrical methods, which can provide 
hints to understand the interfacial magnetoelectricity in 
BiFeOs-based multiferroic heterostructuresJ^i^'2i 

It should be pointed out that the exchange striction 
effect could also induce such a modulation of the bond 
angles, namely the bond angle between a parallel spins' 
pair is different from that between an antiparallel spins' 
pair, e.g. in the multiferroic orthorhombic HoMn03,— 1^ 
which is driven by the non-relativistic exchange (J) be- 
tween spins. This bond modulation can be easily ob- 
served if the lattice is optimized under a collinear G-AFM 
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plus FM configuration: tlie Ru(t)-0-Mn(t) (or Ru(4,)- 
O-Mn(i)) becomes 164.5° while the Ru(t)-0-Mn(i) (or 
Ru(4,)-0-Mn(t)) becomes 168.3°. However, in principle, 
this bond modulation would be switched following the 
switching of FM moments, implying the absence of a 
bias effect although the coercive field might be increased. 
Thus, only the FE polarization can give rise to a robust 
bond modulation with an origin independent of the ex- 
change striction. Of course, the (inverse effect of) ex- 
change striction can affect the underlying FE polariza- 
tion when switching the FM moments, which is also an 
important issue for interfacial magnetoelectricity. 

IV. CONCLUSION 

Summaryzing, in this investigation the ab initio 
method has been used to verify previous proposals that 
were presented to rationalized the exchange bias effect 
in compensated antiferromagnets. The DFT calculation 
discussed here indeed provides qualitative evidence that 
both the Dzyaloshinskii-Moriya interaction and the ferro- 
electric polarization can induce exchange bias effects in 
perovskites heterostructures. Their common ingredient 



is the intrinsic bending of the bonds across the interface, 
which is induced by the oxygen octahedra rotations and 
tiltings. It can be argued that the proposed mechanisms 
originate from "uncompensated exchanges" , namely the 
exchanges {D and J) are no longer uniform but mod- 
ulated following the period of the antiferromagnetic or- 
der. The "uncompensated exchanges" emphasized in this 
work have usually been ignored in the previous models 
for the EB effect that mainly rely on uncompensated mo- 
ments at the AFM interfaces. 
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